Detrended fluctuation analysis on the correlations of complex networks 

under attack and repair strategy 
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We analyze the correlation properties of the Erdos-Renyi random graph (RG) and the Barabasi- 
Albert scale-free network (SF) under the attack and repair strategy with detrended fluctuation 
analysis (DFA). The maximum degree k max , representing the local property of the system, shows 
similar scaling behaviors for random graphs and scale-free networks. The fluctuations are quite 
random at short time scales but display strong anticorrelation at longer time scales under the same 
system size N and different repair probability p re . The average degree (k), revealing the statistical 
property of the system, exhibits completely different scaling behaviors for random graphs and scale- 
free networks. Random graphs display long-range power-law correlations. Scale-free networks are 
uncorrelated at short time scales; while anticorrelated at longer time scales and the anticorrelation 
becoming stronger with the increase of p re . 
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I. INTRODUCTION 

Complex networks are an essential part of modern so- 
ciety. Many social, biological, transportation and com- 
munication systems can be cast into the form of complex 
networks, a set of nodes joined together by links indi- 
cating interactions, yj, H H 3 Recently enormous inter- 
est has been devoted to the study of complex networks. 
Among various existing models of complex networks,two 
different model topologies, the Erdos-Renyi model |5j of 
the random graph (RG) and the Barabasi- Albert model 
of the scale-free network (SF), have been widely stud- 
ied. 0,110 

The random graph is constructed starting from an ini- 
tial condition of N nodes and no edges. Then edges are 
added between pairs of randomly selected nodes with 
connection probability p rg . The number of edges con- 
nected to any particular node is called the degree k of 
that node. The average degree of the random graph is 
(k) = Np rg . The scale-free network puts the emphasis on 
the network dynamics and is constructed with the algo- 
rithm of growth and preferential attachment. According 
to the scale-free model, the network grows over time by 
the addition of new nodes and links. A node newly added 
to the network randomly selects m other nodes to estab- 
lish new links, with a selection probability that increases 
with the number of links of the selected node. One of the 
most relevant is given by power-law degree distribution 
P(k), defined as the probability that a randomly chosen 
node has degree k, P(k) ~ A: -7 . 

Recently the security of complex networks to the ran- 
dom failures or intentional attacks has attracted a great 
deal of attention. [3 [H El El El [II The random 
failure is simulated as the deletion of network nodes or 
links chosen at random, while intentional attack as the 
targeted removal of a specific class of nodes or links. In 
order to protect existing networks and design robust net- 
works, we introduced the attack and repair strategy in 



Ref. [16], which will be described briefly in the following 
section. In this paper, we will study the correlations for 
random graphs and scale-free networks under the attack 
and repair str ateg y with the detrended fluctuation analy- 
sis (DFA) [TtI ITsf and investigate the intrinsic dynamics 
in both networks. The advantages of DFA over many 
methods are that it permits the detection of long-range 
correlations embedded in seemingly nonstationary time 
series, and also avoids the spurious detection of apparent 
long-range correlations that are an artifact of nonstation- 

The paper is organized as follows. In Section II we re- 
view the algorithm of the DFA method. In Section III, wc 
apply the DFA method to investigate the correlations of 
the maximum degree k max and the average degree (k) for 
random graphs and scale-free networks. In last section 
we present our brief conclusions. 



II. DFA 

Detrended fluctuation analysis (DFA) is a well- 
established method for determining the scaling behavior 
of noisy data in the presence of trends without knowing 
their origin and shape. Briefly, the DFA method involves 
the following steps. 



(1) We consider a noisy time series u(i) (i — 1, ..,N) 
and N is the length of the series. We integrate the 
time series u(i) and subtract the mean (u), 
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(2) We divide the integrated signal y(k) into non- 
overlapping boxes of equal size n (scale of analysis). 

(3) In each box of size n, we fit the integrated time 
series y(k) by using a polynomial function, yfit{k), 
which is called the local trend. For ordcr-^ DFA 
(DFA-1 if I = 1, DFA-2 if I = 2 etc.), border poly- 
nomial function is applied for the fitting. 

(4) We detrend the integrated time series y(k) by sub- 
tracting the local trend yjit(k) in each box, and we 
calculate the detrended fluctuation function 
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Y(k)=y(k)-y fit (k). 



(3) 



(5) For a given box size n, we calculate the root-mean- 
square (rms) fluctuation 



F(n) = 



1 N 



(4) 



fc=l 



and repeat the above computation for different box 
sizes n (different scales) to provide a relationship 
between F{n) and n. 

A power-law relation F(n) ~ n a between F(n) and 
the box size n indicates the presence of scaling. The 
parameter a, called the scaling exponent or correlation 
exponent, represents the correlation properties of the sig- 
nal. If a = 0.5, there is no correlation and the signal is an 
uncorrelated signal (white noise); if a < 0.5, the signal is 
anticorrelated; if a > 0.5, there are positive correlations 
in the signal. 



FIG. 1: The maximum degree fc„ 
(b) scale-free network under N = 



ix for (a) random graph and 
1000 and p re = 0.01. 



• Repair: Reconnect this node with the other nodes 
in the network with repair probability p re . 

Then, the evolution comes into the next Monte Carlo 
time step. At each Monte Carlo time step s, we record the 
maximum degree k max {s) and the average degree (k)(s) 
of the system. The reason we choose k max and (k) to 
analyze is that they represent the local and statistical 
features of the system. Since the time step should be 
large enough to enable the system to reach stationary 
state, it is chosen to be 1.5 million in our simulation. 
To reduce the effect of fluctuation on calculated results, 
for every system with size N, the calculated results k max 
and (k) are averaged over 10 independent runs for each 
network realization. 



A. DFA on the maximum degree k n 



III. DFA ON CORRELATIONS OF RG AND SF 

First we will give a brief description about the attack 
and repair strategy. We start by constructing a net- 
work according to the algorithms of Erdos-Renyi random 
graph model (RG) or Barabasi-Albert scale-free model 
(SF). The dynamics of the attack and repair strategy is 
defined in terms of the following two operations: 

• Attack: Find a node with the maximum degree 
kmax and remove all its links. (If several nodes 
happen to have the same highest degree of connec- 
tion, we randomly choose one of them.) 



According to the attack and repair strategy, it can be 
easily seen that the network has a tendency to decrease 
the maximum number of connections among the nodes 
at a long time scale, because the nodes with highest con- 
nections will be damaged and reconnected randomly. 

In Fig. 1, we plot the maximum degree k max versus 
time step s for RG and SF with N = 1000 nodes and re- 
pair probability p re = 0.01. It should be illustrated here 
that, in order to reduce the number of parameters in the 
model, we set the repair probability p re the same as the 
connection probability p rg in RG. Fig. 1 shows that k max 
of RG decreases slowly and fluctuates widely, while k max 
of SF decreases very steeply at the beginning and then 
very slowly with small fluctuations. The different behav- 
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FIG. 2: Scaling behavior of the root-mean-square (rms) fluc- 
tuation function F(n) as a function of the scale n, for series 
kmax of random graphs (RG) and scale-free networks (SF) 
under N = 1000 and p re = 0.01, 0.02, 0.03, respectively. The 
dashed line represents the fit of the slope. 



iors of k max for RG and SF are rooted in their different 
topological structures. RG has homogeneous Poissonian 
degree distribution: all nodes have approximately the 
same number of links. The removal of each node causes 
the same amount of damage. SF is an extremely inho- 
mogeneous network: because the power-law distribution 
implies that the minority of nodes have highly connected 
nodes. The removal of these nodes drastically decreases 
the degrees of the network. 

In Fig. 2 we present results from the DFA-3 analysis for 
the maximum degree k max of RG and SF. The system 
size is N — 1000 nodes and the repair probability p re are 
0.01, 0.02, 0.03, respectively. In this double logarithm 
graph, the root-mean-squarc (rms) fluctuation function 
F(n) as a function of the scale n, exhibit two distinct 
scaling regimes with the crossover at about n = 10 4 , 

F(n)~i n " S n<1 ° 4; (5) 
* W (n* n>10 4 . (5) 

For n < 10 4 , the scaling exponent ot\ is about 0.51; while 
for n > 10 4 , the exponent a 2 is about 0.05. Our DFA re- 
sults suggest that the fluctuation is quite random (close 
to white noise) for short time scales. As the time scale 
becomes larger, the fluctuation displays strong anticorre- 
lation, that is, a large value (compared to the average) of 
k max is almost always followed by a small value. With the 
increase of the repair probability p rei we find that there 
is a vertical shift of the rms fluctuation function F(n) to 
larger values and the shift between p re = 0.01 and 0.02 



FIG. 3: The average degree (k) for (a) random graph and (b) 
scale-free network under N = 1000 and p re = 0.01. 



is very big. In addition, the similar shape between F(n) 
and n is observed for both RG and SF under different 
p re , which illustrates that the time series of k max for RG 
and SF have the similar correlation properties. 

B. DFA on the average degree (k) 

The average degree (k) is a statistical topological prop- 
erty of complex networks. The time series of (k) reflect 
the general changes of RG and SF under the evolution of 
attack and repair. We find that (k) of both RG and SF 
decay exponentially with small fluctuations as shown in 
Fig. 3, under the system size N = 1000 and the repair 
probability p re = 0.01. The inset panel of Fig. 3 gives 
the series of (k) in linear coordinates, showing that (k) 
decreases very steeply at first and then becomes steadily. 
In fact, the curve of (k), to a large degree, depends on 
the repair probability p re . The bigger the repair prob- 
ability, the smoother the curve will be. Here we choose 
p re = 0.01 to analyze, considering to have a comparison 
with the results of k max . 

We next investigate the correlation properties of series 
(k). Fig. 4(a) shows the results of applying DFA-3 to the 
series (k) of RG under N = 1000 and p re =0.01, 0.02, 
0.03, respectively. In the log-log coordinates, the rms 
fluctuation function F(n) and the scale n has a power- 
law dependence with scaling exponent a = 0.64, which 
indicates the presence of the positive correlations and the 
weak long-range (long-memory) power-law correlations. 
The results suggest that a value of (k) at a given time 
depends on values of long past in a power-law fashion. 
We also note that F(n) shifts downwards slightly with 
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law correlations of RG are probably due to the strong re- 
striction, denned as the connection probability equals the 
repair probability, p rg = p re . For scale- free networks, the 
repair probability p re has a great effect on the power-law 
degree distribution of SF in small ranges. The random- 
ness are dominated at small time scales, characterized 
by uncorrelated series. While at large time scales, the 
power-law degree distributions are dominated, character- 
ized by the anticorrelation. 

Further, Fig. 2 and Fig. 4 represent the DFA results 
of series k max and (k). We note that, with the increase of 
p re , the rms fluctuation function F(n) increases in k max 
but decreases in (k). It is the consequence that more 
random connections are added with increasing p re . Ac- 
cording to the local property of k max , the fluctuation of 
the system increases with added randomness. While ac- 
cording to the statistical property of (k), the fluctuation 
of the system decreases since the added connections bal- 
anced the degrees of the whole network. 



FIG. 4: Scaling behavior of the root-mean-square (rms) fluc- 
tuation function F(n) as a function of the scale n, for series (k) 
of (a) random graphs (RG) and (b) scale-free networks (SF), 
under N = 1000 and p re = 0.01, 0.02, 0.03, respectively. The 
dashed line for RG and the solid line for SF represent the fit 
of the slopes. 



the increase of p re , which is oppositive to that of k max 
in Fig. 2. Besides, the big shift between p re — 0.01 and 
0.02 vanishes here. 

Fig. 4(b) reports the DFA-3 results of time series (k) 
for SF under N = 1000 and p re =0.01, 0.02, 0.03, corre- 
spondingly. The two-segment power-law relationship is 
observed with the crossover at about n = 10 3 , 

„, s ( n ai , n< 10 3 ; 
F ( n )~(n% n>10< ® 

For short time scales (n < 10 3 ), the rms fluctuation func- 
tion F{n) has the similar scaling exponent with a\ = 0.55 
for different repair probability p re . For longer time scales 
(n > 10 3 ), the exponents a-2 oip re = 0.01, 0.02, 0.03 cor- 
respond to 0.54, 0.36, 0.23. It shows that the rms fluctu- 
ation function F(n) decreases and the series (k) becomes 
strong anticorrclated with increasing p re for longer time 
scales. 

Fig. 4(a) and 4(b) demonstrate that the rms fluctua- 
tion function F(n) of series (k) exhibits different scaling 
behaviors for RG and SF. The positive long-lange power- 



IV. CONCLUSIONS 

In summary, we use DFA method to analyze the cor- 
relation properties of random graphs and scale-free net- 
works under the attack and repair strategy. According 
to the analysis of the maximum degree k max represent- 
ing the local property of the system, we find that both 
random graphs and scale-free networks exhibit the simi- 
lar scaling behaviors, which suggest that the fluctuations 
are quite random at short time scales but strong anti- 
correlated at longer time scales. With the analysis of the 
average degree (k) representing the statistical property of 
the system, we observe that random graphs and scale-free 
networks display completely different scaling behaviors. 
Random graphs show long-range power-law correlations; 
while scale-free networks arc uncorrelated at short time 
scales but display anticorrelations at longer time scales. 
With the increase of the repair probability p re , the ap- 
parent vertical shifts of the rms fluctuation function F(n) 
for k max and (k) arc due to the effects of the increase of 
random connections to the networks. 
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